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ABSTRACT 

We present the results of a computer simulation of the atomic 
structures of large-angle symmetrical tilt grain boundaries (GBs) S5 
(misorientation angles 36.87° and 53.13°), S13 (misorientation an- 
gles 22.62° and 67.38°). The critical strain level EcrU criterion (phe- 
nomenological criterion) of Chisholm and Pennycook is applied to the 
computer simulation data to estimate the thickness of the nonsuper- 
conducting layer hn enveloping the grain bounclctrics. The h.^ is esti- 
mated also by a bond-valence-sum analysis. We propose that the phe- 
nomenological criterion is caused by the change of the bond lengths 
and valence of atoms in the GB structure on the atomic level. The 
macro- and micro- approaches become consistent if the Scrit is greater 
than in earlier papers. It is predicted that the symmetrical tilt GB 
S5 9 = 53.13° should demonstrate a largest critical current across the 
boundary. 

PACS: 74.50. +t; 74.72. -h; 61.72. Mm 

The major problem in the applications of YBa2Cu307_5 (YBCO) crystals 
is their low critical current density Jc caused by poor current transmission at 
the grain boundaries (GBs). The first studies of this phenomenon, starting 
with |l[] , showed that Jc across the GBs decreases drastically with increasing 
misorientations angle 9 between grains. Subsequent studies have indicated 
that certain specific large-angle GBs (near-Sl 90° [010] GBs) do not show 
this sharp decrease (for details see 0). Several mechanisms have been sug- 
gested to explain the GBs influence on superconducting properties of YBCO 
(see for details 0, Q). We will examine the explanations based on the reduc- 
tion of the order parameter caused by strain. Chisholm and Pennycook ||^ 
suggested that cores of dislocations forming low-angle GB are nontransparent 
for the supercurrent because superconductivity is suppressed in the region 
of the crystal lattice where strains achieve some critical value Ecru- Based 
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on this idea and using experimental data they derived the radius of non- 
superconducting region surrounding dislocation as ~ 2.9|b| (there |b| 
is the magnitude of the Burgers vector of the dislocation). On this point the 
authors of left the semi-microscopical approach and used a strain of 1 % as 
the cut off between YBCO being superconducting and non-superconducting. 
This cut-off value was taken from macroscopic consideration: from the fact 
that, as it is noticed in 0, a 1% strain causes YBCO to be tetragonal and so 
nonsuper conductive. This phenomenological criterion allows one to estimate 
the thickness of the non-superconducting layer enveloping the GB if the dis- 
tribution of the strain in the vicinity of it is known. These type of estimates 
were used for low-angle grain boundaries within the dislocation model (see, 
for example, 

A precise physical meaning of the phenomenological criterion of suppres- 
sion of the superconductivity by strain is not well understood. Within the 
approach p the strain shifts the chemical potential of the Josephson junction 
(JJ) barrier at GB, and this may strongly affect the Josephson coupling, as 
long as the barrier is characterized by the proximity to the metal-insulator 
transition. Many properties of the low-angle GB JJ have been successfully 
explained by this model [Q. It has been pointed out in 0, that small changes 
of the local structure, existing in the vicinity of the GB on the microscopic 
level, should induce large variations of the valence of the copper ions, which 
control the transport properties. This may give a plausible microscopic ex- 
planation for the exponential drop in the critical current with grain boundary 
misorientation 0. 

We propose that the phenomenological criterion of suppression of the su- 
perconductivity by strain ^ is caused by the change of the bond lengths, 
valence of atoms, and the number of charge carriers present in the GB struc- 
ture on the atomic level. Therefore we will try to relate directly the bond- 
valence-sum analysis to the phenomenological strain criterion of suppression 
of the superconductivity. We use a very simple model and assume that all the 
atomic bonds in the given point of the strained lattice experience the same 
change of length irrespective of their orientation. Let us now consider the 
strain of the interatomic bonds in the vicinity of the GB and corresponding 
change of valence caused by it. Each bond between atoms i and j is associ- 
ated with the bond valence Sij that can be determined accordingly to by 
the expression: Sij = exp[(ro — rij)/B], where r^j is the bond length, tq and 
B are empirical constants. The oxidation state of atom i (atom's valence) 
can be written as Vi = Y2j ^ij- If there is the average strain e in the region 
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of lattice in question, the bond length changes and can roughly be estimated 
as r^j ~ rjj(l + e). The atom valence also changes to Vf. It is easy to show 
that e, Vi, and Vf are related by the following approximate relationship: 



s^-llnYL (1) 

Now we can estimate Ecru- If "we take according to a valence cor- 
responding to the superconducting state as = 2.2, and a valence corre- 
sponding to the non-superconducting state as V^'' = 1.8, according to |^ take 
B = 0.37, and bond length between atoms Cul and 04 as r^j ^ 1.96A, we 
will get the following approximate estimation Scru ~ 0.038, that is, several 
times greater than the value of the phenomenological criterion derived in 
§ from the macroscopical consideration. Nevertheless these two approaches 
become consistent if we will come back to the semi-microscopical level of 
consideration. The strain field of dislocation can be estimated accordingly to 
T0[ as e ~ ^ where r is a distance from geometrical center of dislocation. 



If we take r = ~ 2.9|b|, we would get Ecru — 0.027. Thus we can deduce 
that the phenomenological criterion for the suppression of superconductivity 
by strain at lattice level in YBCO should be approximately 3% — 4%. The 
approach described above can be used to estimate the nonsuper conducting 
region in the vicinity of any type of the crystal lattice defects whose strain 
field is known. Therefore it is important to know the precise number of the 

At the misorientations where the dislocation cores begin to overlap, the 
dislocation model of GBs becomes inapplicable (the concept of individual 



dislocations in the boundary becomes meaningless [0). Attempts to use 
the dislocation models to describe large-angle regime are not justified 0. 
Therefore there is no analytical method that allows a consistent calculation 
of the strain field for the large-angle GBs and then to estimate analytically 
the thickness of the nonsuperconducting layer at GB. 

In this work we apply the approach of Chisholm and Pennycook for 
the most perfect large-angle GBs with smallest S (S is inverse of the den- 
sity of coincident sites if lattices of neighboring grains are assumed to fill 
all space). Specifically, we get quantitative information about the atomic 
structure of these GBs by performing their computer simulation and deter- 
mining the distribution of the strain in the vicinity of GBs. Then we use the 
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phenomenological criterion to determine the thickness of the nonsupercon- 
ducting region enveloping large-angle GBs. The thickness of this region is 
estimated also by another method: determination of valences of atoms Cul 
at the vicinity of GBs by analogy to p). 

An inter-atomic potential for YBa2Cu307 based on an unscreened rigid- 
ion model has been developed in This potential was used for computer 
simulation of twin boundary (TB) |]l3l and TB and GB [l^ in YBa2Cu307. 
In |jl5|, the empirical interatomic potential model for YBa2Cu307_5 was pro- 



posed. We employ this potential |jT5| in this study as well as the Ewald 
method of calculation. 

We will characterize GBs by the direction of the axis of misorientation, 
the angle of misorientation relatively to this axis and the geometrical plane of 
GB which is chosen coinciding with some simple crystal lattice plane. In this 
article we present results of the computer simulation of the atomic structures 
of the symmetrical large-angle tilt GBs with the misorientation axes [001]: 
S5 (310)1/(310)2, misorientation angle 6 = 36.87°; S5 (210)i/(210)2 9 = 
53.13°; S13 (510)i/(510)2 9 = 22.62°; S13 (320)i/(320)2 9 = 67.38°; S29 
(520)1/(520)2 9 = 43.60°. Boundary planes common for the two grains are 
given in parenthesis. The indices 1 and 2 refer to the two neighboring crystals 
of bicrystal. In general, the plane which is shared by the two crystals can 
be described by specifying the planes in each of the two crystals which is 
parallel with the common boundary plane |jl6[. The GBs S5 and S13 were 
chosen as the most perfect tilt GBs that cover wide range of misorientation 



angles and are experimentally observable [|17|, |16|. The S29 GB is chosen as 
an example of GB close to the GBs of common type. 

The configurations of models are following. The XOY-plane coincides 
with the basal ab-plane of the lattice, the XOZ-plane coincides with the 
geometrical plane of GBs. The X-axis lies at the GBs geometrical planes. 
The Y-axis is perpendicular to the GBs geometrical planes. The Z-axis is 
directed along c-axis of crystal lattice and coincides with the misorientation 
axis. The extent of models along X-axis in all cases are equal to one period 
of coincident-site-lattice of corresponding GBs. The total number of ions in 
model bicrystals for GBs S5 ^ = 36.87°, E5 ^ = 53.13°, E13 9 = 67.38° is 
624. For GBs S13 9 = 22.62°, S29 9 = 43.60° it is 650. 

Periodic boundary conditions were applied at all the outer faces of the 
computational cell. The molecular dynamics method was used. The equa- 
tions of motion were solved using the velocity form of the Verlet algorithm 
with the starting time step not greater than 2, 5 ■ 10~^^ s. In each of cases 
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of GBs simulation, we started with a sharp geometric initial configuration of 
the boundary. It took us 224 time steps to achieve relaxed configuration in 
the case of GB S5 ^ = 36.87°, 122 in the case of GB S5 ^ = 53.13°, 2471 
time steps in the case of GB S13 9 = 22.62°, 4331 time steps in the case 
of E13 e = 67.38°. In the case of GB S29 6 = 43.60° the system was far 
from equilibrium after 3691 time steps and there are serious doubts that the 
equilibrium is achievable with a sharp geometric initial configuration of the 
E29 GB. This means that, for such a boundary, the equilibrium could be 
reached only if a concentration of vacancies exists at the GB. 

Relaxed atomic structure of GB E5, 9 = 53.13° is represented in the 
Fig. 1, GB S13, 9 = 22.62° in the Fig. 2 (all distances are expressed 
in angstroms). The layer representation of the YBCO lattice cell is 
used. All layers are parallel to the basal plane. The plane consisted of Y 
atoms is a refiection plane. Thus it is easy to imagine the atomic structure 
of rest lattice along Z-axis. Comparison of the initial and final (relaxed) 
configurations allowed to determine the atomic displacement vectors for the 
all atomic sites. They are shown as arrows in the figures. It was found 
that for all GBs, the z-components of the atomic diplacement vectors are 
significantly smaller than x- and y-components. The maximal displacements 
are in the atomic plane next to the geometric plane of the GB. Especially large 
displacements have oxygen atoms of types 04 and 03. The magnitude of 
displacements drastically drops with distance from the GB geometric plane. 

In order to facilitate analysis of spatial displacement distribution the aver- 
age atomic displacement vector u was found in the each atomic plane parallel 
to the geometric plane of GB. The components of this vector were found by 
averaging of corresponding components of all the atoms at this atomic plane. 
For GB S5 = 36.87°, the magnitude of the average displacement vector 
u = |u| at the plane of the GB is 1.54 ■ 10~^ A. At the adjacent to the GB 
plane, u = 5.53 ■ 10~^ A, at the distance from GB y = 4.23 A it drops by the 
order of magnitude, at the distance from GB y = 6.64 A it drops by the two 
orders of magnitude and so on. The same character of displacement distribu- 
tion is occured for GBs S5 = 53.13° and S13. The smallest displacements 
are in the vicinity of GB E5 ^ = 53.13°, the greatest in the vicinity of GB 
S13 9 = 67.38 °. In the case of GB S29 the displacement of atoms 03 and 04 
from sites adjacent to the geometrical plane of GB are extremely large 3 A 
and there is no evidence of their slowing down - configuration far from relaxed 
after 3691 time steps calculation. It seems that in relaxed configuration we 
would expect depletion of oxygen atoms at the GB. 
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Our results agree with experimental observations of the atomic structure 
GB E5 ^ = 36.87° that the symmetric S5 ^ = 36.87° GB demonstrates 
no detectable depletion zone and the lack of disordered region. The displace- 
ments of oxygen atoms outward of GBs and drastic increase of relaxation 
time for S29 allows to suggest that the GBs with large E could exist only 
if the oxygen depletion zone is formed in the vicinity of GBs. A presence of 



these zones, as it have been experimentally shown in [19, 20|, may essentially 



influence the superconducting properties of the large-angle GBs. 

Now we will try to relate atomic structure and transport properties of 
GBs estimating hn from computer simulation data. The results of the bond- 
valence-sum analysis applied to the relaxed atomic structures of the GBs ob- 
tained by the computer simulation are represented in the Table. The distance 
from the GB plane where bond- valence sum becomes close to the same in the 
"bulk" is denoted as r„ {hn = 2r„). Then we will apply strain analysis |n . 
The effect does not significantly depends on the type and sign of strain 0. 
Therefore to facilitate analysis only the absolute value of the average atomic 
displacement |u| in each atomic plane parallel to GB were used to construct 
approximate functions |u(r)| (r is the distance from the geometric plane of 
the GB). A measure of the lattice strain e was estimated as the derivative of 
the approximate function with respect to r. The values of hn corresponding 
^crit = 3% are represented in the Table. Analysis of data shows for all cases 
that the smallest /i„ has GB S5 = 53.13°. The GB S5 = 36.87° has 
slightly greater /i„ though its misorientation angle is smaller. The hn of other 
considered GBs increases with the increasing of misorientation angles. The 
width of non-superconducting zone hn can be used to estimate the expected 
tunneling current Jc that exponentially decreases with hn- One can expect 
that the more perfect special tilt GB S5 = 53.13° characterized by the 
smallest length of the periodicity along GB should demonstrate the larger 
Jc across boundary. We can suppose that this GB has the best transport 
properties in comparison with other symmetrical tilt [001] GBs. The GB E5 
6 = 36.87° should have smaller Jc though its misorientation angle is smaller. 
The two symmetrical GBs, which are complementary to each other, are dif- 
ferent on an atomic scale and therefore, as it was pointed out in [1^, the 
electrical measurements need to be carried out on the whole range from ° 
to 90°. Jc for the other analyzed GBs should decrease with the increasing 
of misorientation angle in agreement with the experiment. The results show 
that suggestion of suppression of superconductivity at GBs by strain leads 
to very sharp difference in Jc through special GBs with small E. We sug- 
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gest that this prediction could be tested experimentally (these results are 
relevant only for unfaceted GBs) and could shed light on applicability of the 
mechanism of suppression superconductivity by strain at the GBs in high-Tc 
sup er conduct ors . 
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Table 1: Thicknesses of nonsuperconducting layers enveloping GBs with mis- 
orientation axis [001] in YBa2Cu307_5. 
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Figure 1: Atomic structure of GB S5, 9 = 53.63°. 

Y - +, Ba - X, Cul - ^, Cu2 - oP- ■, 02 - 0, 03 - •, 04 - A. 

The layer ^1 (atoms Cul and 04)is shown in Fig. 1(b). The layers ^^2 

(atoms 01), #3 (atoms Ba) are shown in Fig. 1(c). The layers #4 (atoms 

Cu2), #5 (atoms 02), #6 (atoms Y) are shown in Fig. 1(d). The projections 

of layers #1 — #6 are shown in Fig 1(a). 
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(a) 
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Figure 2: Atomic structure of GB E13, 6 = 22.62°. 
The designations are the same as in Fig. 1. 
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